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Abstract 

We diagonalize the transfer matrix of a solvable vertex model constructed by com- 
bining the vector representation of U q [Sl(n\m)] and its dual by means of the quantum 
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to derive the transfer matrix eigenvalues and the respective Bethe ansatz equations for 
general grading choices. 

PACS numbers: 05.50+q, 02.30.IK 
Keywords: Algebraic Bethe Ansatz, Lattice Models 
December 2005 



1 Introduction 



Nowadays it is well known that the concept of integrability based on the Yang-Baxter alge- 
bra pQ accommodates in its algebraic structure variables with commuting and anti-commuting 
rules of permutation [2 . This framework provided us a systematic way of constructing and 
solving two-dimensional integrable vertex models possessing both bosonic and fermionic edge 
configurations. The respective Boltzmann weights are viewed as i?-matrix solutions of the 
graded version of the Yang-Baxter equation based on supergroup symmetries The associ- 
ated transfer matrices can then be diagonalized algebraically by means of the quantum inverse 
scattering method [Bj. 

An important class of the graded vertex system is that directly related to the finite dimen- 
sional representations of the Sl{n\m) superalgebra j3]. The corresponding i?-matrices have 
been given origin to a number of interesting one-dimensional magnets and interacting elec- 
trons systems. To name just some examples, we mention the lattice gas model with n species 
of bosons and m species of fermions [3], the supersymmetric t — J model [HJ U\, integrable 
generalizations of the Hubbard model [HI E] and Heisenberg spin chains JU] • 

Integrable impurities ^Tj can also be introduced into these systems by considering different 
representations of the Sl(n\m) symmetry at some lattice sites. In this manner, an impurity 
t — J model [H] and a family of doped Heisenberg chains [T3j, related to four- dimensional 
representations of S7(2|l), have been proposed and solved by algebraic Bethe ansatz. One way 
of incorporating defects with the same Hilbert space as the bulk is by combining the vector 
and the dual representations of Sl(n\m). This possibility has been explored not only in the 
context of impurity vertex models ^3] but also to model an integrable network system to 
study the corner transfer matrix spectrum ^Hl- However, as far as Bethe ansatz results are 
concerned, such construction has been restricted to the specific case of the rational Sl(2\l) 
vertex model ^3]. This solution has recently been important to understand the universality 
class of a spin chain that alternates between the 3 and 3 representations of Sl(2\l) |17j . 

These considerations indicate that the transfer matrix diagonalization of the vertex model 
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whose weights involve both the vector and dual representation of the extended U q [Sl(n\m)] 
quantum superalgebra, for arbitrary n and m, is worth to be pursued. This model contains 
an extra free parameter q in such way that the classical Sl(n\m) symmetry is recovered at the 
isotropic q — > 1 limit. Here we tackle this problem by means of the algebraic Bethe ansatz such 
that the many possible grading choices are considered in a unified way for arbitrary n and m. 
We recall that the latter feature for general values of the number of bosons and fermions has 
been known to be technically involved in the literature [HI El EI]- The typical difficulty one finds 
is to establish a suitable formulation of the nested Bethe ansatz that takes arbitrary grading 
into account and therefore avoiding a case by case analysis |18j . In this paper this problem is 
circumvented and the transfer matrix eigenvalues as well as the corresponding Bethe ansatz 
equations are exhibited for any of the ^i^yl possible grading choices. 

This paper is organized as follows. In the next section we describe the mixed vertex model 
whose weights alternates between the vector representation of U q [Sl(n\m)} and its dual. The 
corresponding Hamiltonians representing either impurity or alternating spin chain are explicitly 
derived. In section E] we describe the essential steps to solve the eigenvalue problem for the 
respective mixed transfer matrix. In section 0] explicit expressions for the eigenvalues and 
Bethe ansatz equations are presented for general grading. In Appendix A we summarize the 
coefficients of the Hamiltonians given in section El 

2 The vertex model 

The system we are going to consider is determined by the fundamental i?-matrix associated 
with the vector representation of the U q [Sl(n\m)] superalgebra ^H]- We recall here that this 
i?-matrix is closely related to the Boltzmann weights of the so-called Perk-Schultz vertex model 
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and it is given by 



n+m n+m 

a=l ct,/3=l 
n+m n+m 



+ d(A) £ g S ® g $ + c 2 (A) £ 4 a 2 ® 4% C 1 ) 

a</3 a>/3 

where the Boltzmann weights a a (A), 6(A) and Cj(A) are determined by 

g2(l-p a ) _ q 2p ae 2\ 

a a (X) = _ 2A , a = l,...,n + m, (2a) 

6(A) = qi l~ e l X x \ (2b) 
Ci (A) = [q g2 j e2A , , = 1,2 (2c) 

where G C" +m are the standard (n + m) x (n + m) Weyl matrices. The symbols q and 
A denote the quantum group deformation parameter and the spectral parameter, respectively. 
The Grassmann parities p a distinguish the bosonic p a = and fermionic p a = 1 degrees of 
freedom. We stress that in this paper the grading ordering of the many possible permutations 
of pi . . . p n+m will be considered arbitrary. 

The graded Yang-Baxter i2J satisfied by the i?-matrix ([TJ is the basic object to build 
the transfer matrices of solvable vertex models mixing between different representations of 
U q [Sl(n\m)]. This is an associative algebra generated by the elements of the Lax operators 
C A j(\), namely 

R(\ - ii)C Aj {\) & C Aj {y) = C Aj (fx) J> £ Aj {\)R{\ - n), (3) 

s 

where the supertensor product A <g> B between two matrices with elements A a p and is 

n+m 

viewed as (A I B) X2 = £ (-l) p ^+^A aj B ps e^ <g> ef 8 . 

o/37<5=l 

The operators C A j(X) are matrices on the auxiliary space A = C n+m whose elements are 
operators on the quantum space Vj G U q [Sl(n\m)]. The simplest Lax operator occurs when Vj 
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is chosen to be the space of the vector representation of U q [Sl(n\m)]. In this case the algebra 
(J3J) coincides with the Yang-Baxter equation for R a b(X) provided we set 

4 b \\)=P ab R ab (\), (4) 

n+m 

where P a b is the graded permutation operator P a b = (— l) PaPl3 e^ ® e^ a . 

a,/3=l 

The other possibility we are interested in is the situation in which Vj is taken as the dual 
of the vector representation. The general theory concerning the construction of universal 
Lax operators for the U q [Sl{n\m)] superalgebra have already been discussed in ref. [21]. To 
obtain from that suitable expressions for a posteriori Bethe ansatz analysis requires, however, 
a certain amount of additional work. Omitting here these steps we find that the Lax operator 
intertwining between the vector and conjugated representation of U q [Sl(n\m)} in the Weyl basis 
is 

n+m n+m 



3 W 

-A3 

o=l a,0=l 



n+m 



+ ci(a) (-^- 25a < a) , ® €1 + m e (-w^- 1 ^ ® ^, ( 5 ) 



a,3=l a,/3=l 
Q</3 a>/3 



where the new Boltzmann weights a a (X), 6(A) and q(A) are given by 



a Q (A) = — , a = l,...,n + m, (6a) 



1 -e 2 
1 {q 2 - e 2 A ^ 
g 1 — e 2 



ka) = ; 2A ; , (6b) 



Ci{\) = ^^e 2 ^ 1 ), i = l,2. (6c) 



Having found two different Lax operators satisfying the Yang-Baxter algebra with the 

ing framework. As usual the corresponding row-to-row transfer matrix T^ Ll,L2 ^(A) is written 
as the supertrace over the auxiliary space of an ordered operator called monodromy matrix 



same i?-matrix, an integrable vertex model that combines L\ operators of type £^](A) and L2 
operators of type C A AX) can be directly build up by means of the quantum inverse scatter- 



Tj Ll ' L2) (A), namely 



T (Ll ' La) (A) = str^ 



n+m 



i=l 



(A) 



where the monodromy matrix is given by 



(7) 



^4 Ll ' i2) (A) — C-AL 1 +L 2 W£-AL 1 +L 2 -l(^) ■ ■ ■ Aa2(A)Aai(A). 



The new operator C A j(\) is now defined as follow 



C Aj (X) 



^S](A), 



if 3 e {&,..., 



£S](A + p), if./f U, 0,}. 



(9) 



where the partition {£i . . . ^} denotes a set of integers indices assuming values on the interval 
1 < £j < L = Li + L 2 . The inhomogeneity p can be arbitrarily chosen and it does not spoil 
basic properties such as the symmetry and locality of the interaction. 

Next we turn our attention to interesting types of quantum spin chains that can be derived 
from the transfer matrix 0. Though the integrability does not depend on how we choose the 
partition . . . £z, x } the construction of local conserved charges commuting with T^ Ll ' i2 )(A) 
does. One interesting situation occurs when we have only one operator of type C\ L (X) sitting 
at the end of the chain. This plays the role of an impurity on a chain of L — 1 C^(X) Boltzmann 
weights P3] . The respective impurity Hamiltonian is obtained as the logarithmic derivative 
of the transfer matrix ((7|) at the regular point A = 0, and reads 



L-2 



= E 4+i (°) 



-i 



A ( ? + i(o)+ 41» cg-iM 



-1 



1(2) 



i=l 



(10) 



To obtain explicit expressions for this operator one needs to known the form of the in- 

(2) 

verse of £ A j(\). Since this operator combines different representations, its inverse does not 
follows directly from unitarity. After some cumbersome algebraic manipulations we obtain 
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that 



C j\j{\) is given by 



4? (A) 



n+m 



n+m 



a=l 



a,/3=l 



n+m 



n+m 



(11) 



+ c 1 (A)E(-l) Pa C a yi>e2 + c 2 (A)E 

a>/3 



-e a/9 <x> e a/3 , 



a, 0=1 

a</3 



where the functions a a (A), 6(A) and c,(A) are 



6(A) 
a (A) 



2 a^) (l-^ A )(l- e2 V (n " m ^ 2+2pa) ) 

^ (g2 _ e 2X\(\ — e 2\q2(n-m-l)\ 

g(l - e 2A ) 
(g 2 - e 2A ) ' 

(g 2 -l)(l-e 2 y(— g 2A(t _ 1} 

(g 2 -e 2A )(l - e 2A g 2(n- m -l)) e 



a = 1, 



, n + m, 



z = l,2, 



and 



^ = g -2(i-i+5i >1 ) g 2(P i +P^)pj ? -4p &) 



(12a) 
(12b) 
(12c) 

(13) 



k=i 



for i < j, i, j — 1, . . . , n + m. 

Substituting Eq. (fTT|) into Eq. ffTUj) and after several simplifications we find that US 1 ' can be 
rewritten as 

^ (1) = E^i + ^ + ^w ( 14 ) 



i=l 



where h^- +1 , hfj +1 and /ij^^j 1 are defined as follows 



h 



(i) 

(2) 



-Ri,i+l(0), 
n+m 



(15) 



n+m 



n+m 



KU = E ® + E ->fo« ® + E e «^$ ® ^ +1) (i6) 



a,/3 = l 



a,/3=l 
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n+rn 



a=l 
n+m 

+ £ 6(p)a a (0)5(p)(-l)^e£r 1) ® 4? ® 32 

a,/3=l 

n+m n+m 

+ E ® g 2 ® $ + E ^Sr 1 ' ® $ ® 42 

a,/3=l a,/3=l 
n+m n+m 

+ E ® ® $ + E wSr 1} ® «S ® 45 m 

a,/3=l a,/9,7=l 
n+m n+m 

ct,/3,7=l a,/3,7=l 
n+m n+m 

a,/3=l a,/3,7=l 

such that the coefficients fi^g, and S a)( g i7 are exhibited in Appendix A. The diagonal 

matrix M a/ g has the following form 

L-3 

n+m (Pa+Pff) / ] Pv k 

Maj,= E ^ ^l®41®---®4t 3 l-3- ( 18 ) 

I/l,...,l/ L _ 3 = l 

Yet another relevant case occurs when the number of operators Cj^AX) and Cj^AX) are 
equally distributed Li = L2 = L/2 in an alternating way in the monodromy matrix [221 Ell- 
in this situation the respective Hamiltonian is a bit more involved and is given by 

n {2) = E [^Up)\ 1 ^Up)+ E [4i\>)]~ X [^ 

(19) 



i=l i=2 
odd i even i 



Its explicit expression in terms of the Weyl basis, after some algebraic simplifications, 
involves many of the terms entering in the previous Hamiltonian. More precisely, the final 
result is 



z=l i=2 

odd % even i 
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— (ON 

where the new bulk term h i _ 1 i i+l is given by 

n+m 

E a Q (p)a a (0)a Q (p)(-l)^C- 1) ® e£ ® e 
0=1 

n+m 

E &(pK(0)&(p)(-l)^e 



'^-1,1,1+1 



(<+i) 

aa 



a,/3=l 
n+m 



n+m 



(i+l) 
aa 



a,/3=l 
n+m 



a,/3=l 
n+m 



+ il a,/3 e a/3 W e af3 ^ e /3/3 + ^ a Al e - 



a/3 



e a/3 ^ ^ry" 1 



(21) 



a,/3=l 
n+m 



a,/3,7=l 
n+m 



a/3 ^ e 77 ® e /3a 



a,/3,7=l a,/3,7=l 
n+m 



(i+l) 



(Po+p 7 )j) 7 5(*-l) <> -(i+l) 



-a/3 



e 7^ «> e 7 a 



a,/3=l 
n+m 



+ E r2 



(1) 



"a/3 



a,/3,7=l 



From the form of |^C^(A + p) one clearly sees that the above expressions for T^W an d 
7i( 2 ) for p = make sense only when n — m ^ 1. In fact, for generic values of g we observe 



det 



1 — e 2A )£^(A) | a=o vanishes for n — m = 1. It should be remarked, however, that such 



singularity is easily avoided in the rational limit q — > 1. 



3 Algebraic Bethe ansatz 

The purpose of this section is to describe the essential tools entering in the diagonalization 
problem of the transfer matrix (J7J), namely 

T^)(A)|0>=A(A)|0>, (22) 

by means of the nested Bethe ansatz approach 23j. 

An important requirement in this method is the existence of an appropriate reference state 
|0) with the following properties. The action of the monodromy operator on this state gives us 
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a triangular matrix as well as it is an exact eigenvector of the diagonal elements of T^ Ll ' L2 \\). 
The structure of the ^-operators of previous section suggests that such state can be built up 
by the tensor product 

(23) 



i°>= n ®io (i) >,® n ®i° (2 v 



where the local pseudovacuums (i = 1,2) acting on the j-ih site of the L\ + L 2 chain 

are the following vectors with length (n + m) 



|0«) 



, |0( 2 )> 



V / 



o 



(24) 



V 1 / 

(n+m) \ / (n+m) 



It is not difficult to see that the action of each operator C^(X) on its respective local state 
. gives as result an upper triangular matrix, namely 



4](A) |0 (1) > 3 



V 



and 



A (A) # 

/ 2 (A) 






( <7i(A) 
<? 2 (A) 



# 




# 




\ 



/n+m-l(A) 

/n+m (A) y 



|o«y 



(25) 



(n+m)x(n+m) 








# 
# 



\ 



V 








9 n+m— I (A) # 

#n+m(A) y 



|0 (2) >, 



, (26) 



(n+m)x(n+m) 



where the symbol # stands for non-null values. The functions /i(A) and gi(\) are given by 

{ai(A), % — 1, 6(A + p), i = 1, . . . , n + m — 1 

, 9i(\) = \ (27) 
6(A), i — 2,...,n + m a n+m (A + p), i = n + m. 
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The next step is to write a suitable ansatz for the monodromy matrix T^ 1,L2 \\) in the 
auxiliary space A. The above triangular properties suggest us to seek for the structure used 
in the nested Bethe ansatz diagonalization of SU(n) vertex models 



A(X) 
Ci(A) 



B 1 {X) 
Di,i(X) 



,n+m— 1 (A) 



\ 



(28) 



l Cn+m— i(A) -D n + m _i j i(A) • ■ ■ D n j rrn _i n j rm _i 

(A) / 

\ / (n+in)x(n+m) 

and in terms of this representation the eigenvalue problem (J22)) becomes 



n+m— 1 



-1)*M(A) + J2 (-l) ft+1 A,(A) 



|0) = A (A) 



(29) 



An immediate consequence of properties (J25j) and (J26|) is that one can easily derive how the 
elements of the monodromy matrix acts on the reference state |0). The action of the diagonal 
fields turns out to be 



A(X) |0> 
Ai(A) |0) 



Uiwr [giW} L2 io> 



[f i+ x(X)r [g* +1 (\)] L2 |0) , % = 1, . . . , n + m - 1, 



(30) 
(31) 



while that of the non-diagonal operators are 



Bi(X) |0) = #, D i;j (X) |0) = 0, i^j,j^n + m- 1, 
Ci(A)|0) = 0, A,n+m-i(A) |0) = #, + 



(32) 
(33) 



This reveals us that B^X) are creation fields with respect to the reference state |0). The next 
step in the algebraic Bethe ansatz method is to search for other transfer matrix eigenvectors 
as linear combinations of products of such creation fields acting on |0), namely 



B ai (X?)...B a ^(X^ a ^\0) 



(34) 



where sum over repeated indices = 1, . . . , n + m — 1 is assumed and J^ 01 ^- 01 ^ are coefficients 
of the linear combination. 
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In order to solve the problem (|29|) one has to carry on the fields A(X) and Da(X) over the 
operators _Bj(A^ ) until they reach the reference state |0). This is done with the help of an 
appropriate set of commutation rules coming from the Yang-Baxter (J3J) for T^ 1,Ij2 \X). The 
most useful relations are 

A(A)£» = a ^Zx) B ^ )A{ - X) ~ C ^-xj Bj{X)A{fl){ - 1)P1 ' (35a) 

D %j (X)BM = ^±— Bl ^)D iq {x)f^{x - 

- ~ /A . ) ^-(A)Afc(/x)(-l) plp<+1+Pj+l(pi+Pi+l) , (35b) 

0(A - \L) 

Bi{X)B^) = *_ B l { l i)B q {X)f W(A - ^(-1)^+^+0+^ (35c) 
where r^(A)^t are the elements of the i7 g [SZ(n — 1 +j?i|m — pi)] i?-matrix {IJ). Here we remark 

n+m— 1 

that such elements are defined by r^(A) = ^^Wa^cry ® anc ^ the new parities pa 

a[3-y8=l 

defining the auxiliary matrix r^(A) are therefore = for a = 1, . . . , n + m — 1. 

By iterating the fields A(X) and Du(X) over the multiparticle state (jSljl we generate states 
proportional to |</>) as well as those that are not usually denominated unwanted terms. The 
former terms contribute to the eigenvalue A (A) and are obtained by keeping only the first 



term of the commutation rules ()35a|) and (|35b[) . The unwanted terms, however, appear as 
a consequence of the second part of the relations ()35aj) and (|35b|) when the variables A^ 
are exchanged with the spectral parameter A, which can be collected in closed form thanks 
to Eq. (|35c|) . Since this procedure is by now standard in the literature we just list the main 
results. The corresponding eigenvalues are given by 

™i „ t \ (!) _ \ \ mi 

A(A) = {-IT [/i(A)] Ll MA)] La II zkk * + II 77\ ,d)^ A(1) ( A ' ( 36 ) 

j= i b(X)> -X) J=1 b{X-X)>) 

provided that the rapidities {A^ } satisfy the following Bethe ansatz equations 



(37) 
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for i = l,..., fh\. 

The auxiliary function A^(A, {A^}) turns out to be the eigenvalues of an inhomogeneous 
transfer matrix, namely 

r (1) (A, {A^})^:;:;^---- |0) = A«(A, {aW})^"^ |0) , (38) 

whose matrix elements are defined as 

t^\\, {x^})^ 1 '"^ 1 = (— i) ft + i ( i+ sriiP,flfc+i)+piErir 1 p« fe +i+(" i i- i )pip>+if( i )( - x _ ji 

r (1) (A - A£t 1 )&: 2 S::^ (1) (A - Affi^^^w- ( 39 ) 

This completes the first step of the Bethe ansatz analysis since we still need to determine 
the eigenvalues AW(A, {A^ }). In the next section we will discuss how this problem can be 
solved in terms of recurrence relations. 



rewritten as T«(A, {X®}) = str^ (1) 7$, (A, {A^}) 



4 Eigenvalues and Bethe ansatz equations 

Here we are concerned with the diagonalization of the auxiliary eigenvalue problem ()38|) whose 
solution needs a second algebraic Bethe ansatz analysis. The transfer matrix (|39|) can be 

(A, {A^}) with A m = C n+m - 1 such that its corre- 
spondent monodromy matrix is given by 

2$ (A, {Ai 1} }) = D AW (X).r% )rhi (X - xgy^^iX - xgU) . ..r% n (X - A?>), (40) 

where r^ 1} .(A) = PaW^j} 1 ) -(A) an d ^^(^(A) is a (n + m — 1) x (n + m — 1) matrix whose 
elements are the operators Dij(X). We also recall the products in (j4*Uj) are taken in the auxiliary 
space A^. 

An important property of this auxiliary monodromy operator is that it satisfies the following 
intertwining relation 

fW(A - Ai)2$(A, {A«}) § 7$(ai, {A?}) = T%{^ {A< 1J }) § 7«(A, {A^f «(A - p), 

(41) 

12 



S1 (1) 

where now ® stands for the supertensor product with parities p a . 

The next step is to exhibit a pseudovacuum state |0) on which the monodromy matrix 
^l(i)(A, {A^ }) acts as an upper triangular matrix. One way to achieve this is by a suitable 



combination of vectors acting on the tensor product space n^=[ i2 ®C™ +m <8> Yljli ®C 
|24j . Once the first part of this tensor product was already determined as being |0), the second 
one should be set up provided the triangularity property is reached. It turns out that such 
appropriate state is 

I A 





~m+m— 1 



10} 



(1) 



mi 



|o> ^JX 1 



i=l 



(42) 



V / 



(n+m— 1) 



where the monodromy matrix elements of T^i } (A, {Ajj. }) proportional to Dij(X) act on the 
first part of |0)^ while the remaining operators act on the second part of the tensor product 



In fact, by writing (A, {A^ }) in the auxiliary space as 



^gi(A,{Ai 1) }) 



( A« (A) B[ 1] (X) 
cf(A) D«(A) 



\ ^n+m-zi^) ^n+m-2,l(A) ' 

one can derive the following relations 

mi 



D 



(i) 

n+m— 2, n+m— 2 



(A) /. 

/ (n 



(n+m— 1) X (n+m— 1) 



(43) 



(44) 



i=i 



^(A^A^WIO)^ = [/ J + 2 (A)] Ll b+2(A)] L2 nKA-A5 1) )|0) ( 
i = 1, . . . , n + m — 2 



ln\(!) 



(45) 
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^(MA^DIO)^ = 0, i^j,j^n + m-2, (46) 
A^+^tMA^}) |0> (1) = #, + (47) 

cf^MAj?}) |0) (1) = i = l,n + m-2. (48) 

From now on the basic steps in order to diagonalize the transfer matrix T^(X, {X k }) 

become similar to those described in the previous section. One introduces a second mul- 

tiparticle state = B$(\f\ {A^}) . . . B2l 2 (A^, {X^})^^ 2 '"^ 1 |0) (1) parameterized 

by a new set of inhomogeneities x[ . ..A^. Thanks to the Yang-Baxter algebra (f4*T]) sat- 
r (D 



isfied by (A, {A^ }) the structure of the commutation rules between the new diagonal 



fields AW(A,{a£ 1} }) and D$\\, {A^}) are exactly the same as that presented in section El 

At a general level I each operator 0(A) in Eqs. (j35aHH5cj) is substituted by its corresponding 

r i i 

0^(A, {A^}). The auxiliary i?-matrix r^(A) is that of the U q Sl(n — I + ^^Pk\fn — ^^Pk) 

k=l k=l 

vertex model and the parities p^} = p^+i for a = 1, . . . ,n + m — I. As a consequence of that 
the eigenvalues at subsequent steps I and I + 1 are going to satisfy a recurrence relation similar 
to that of Eq. (J3T?|) . More precisely, for a given I > 1 we find that 



A ( ^A,{A?}) = (-ir-[/m(A)] Ll [^ + i(A)] L2 n a m(A-Af) ][ ,+1 ' 



^ +1 a m (Af 1} -A) 
(i+i) 



fhi fhi+i 



Kxr 1 ' - a) 



+ 11^-^) I] 77T \ v+1) . ^ +t> (.\{^,---^}), (49) 

3=1 3=1 »l A - A i ) 

provided that the set of variables {X^ +1 ^} that parameterize the eigenvectors of the inhomoge- 
neous transfer matrix T^(A, {A^, . . . , A^ ; }) satisfy the following Bethe ansatz equations 



/ i+1 (a, (,+i >) 91+1 (as i+ ") 'n : :„ :r n**r ) -*. 



A i ~ A i I TT h (\(l+l) _ )(W)^ ffl '+U A j -\ ) 
A KA," + " - Af ) i\ 1 ' 6(Af " - A?+") 

A( m )(A = A (m }) {A (i ); _ _ _ ^ A (m )}); . = ^ _ ^ (50) 

for I = 1, . . . , n + m — 3. 

Those equations can be iterated starting with I = 1 until we reach the nested Bethe ansatz 
level I = n + m — 3. In the last step / = n + m — 2 one has to deal with a diagonalization of six 
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vertex like model. In this case the final result for the eigenvalue pS n+m 2 \X, {A^ n+m 2 ^}) is 



A (n +m -2) (A; | A (n + m-2) }) 



(51) 



mn+m-2 



{-iy^-i [/ n+m _ 1 (A)] Ll [g n+m -i(X)] L2 n «n +m -i(A-A 



(n+m— 2) >, 
3 1 



•n+m-1 /\ (n+m-1) \\ 

nO'n+m-iy^j — A J 



i=i 



m n+m-2 



(n+m-2)^ 



.7=1 

^n+m— 1 



fo(A ("+m-l) - A) 



t/„ +m (A)] Li K +m (A)] La n 6 ( A " A 



whose Bethe ansatz equations for the last variables {X^ +m ^} are 



n a n+m(A — A^ ■*) 



6(A- A 



(n+m-1) \ 



f n+m-l (A 



(n+m— 1) ^ 



f /\ (n+m-l) x 
Jn+ml^j , 



Li 



5'n+m-l(A 



(n+m— 1) \ 



5'n+m(A 



(n+m— 1) \ 



L-2 



m-n + m-2 



n 



b(X 



(n+m— 1) ->(n+m— 2)\ m ;+1 



n 



n ( V 

"n+mV^i 



(n+m— 1) ^(n+m— 1)\ 



^^ A (™+m-l) _ A (n+m-l)\ 



L J. /\ (n+m— 1) \(n+m— 2)\ J. J. /■> (n+m— 1) \(n+m— 1)\ » / ^ (n+m— 1) \ (n+m— 1)\ ' 

,=1 a n+m-l(A — A ? - J J=1 OVH-m-H^j — Aj J 0(Aj — Aj J 

(52) 



3= 



We have now reached a point in which the results of section |21 and |3] can be put altogether. 
Combining Eq. (J36|) with the iteration of Eq . (J49j) until I = n + m — 3 and finally with the help 
of the expression (J52|) we find that the eigenvalues A(A, {A^. }) are 



A(A; {AW}, . . . , {xf +m - l) }) = (-1)^ [/i(A)] Ll [^i(A)] L2 \\ 



m ai(A; i} -A) 



5 6(A« - A) 



n+m— 1 



(n+m— 1) ^ 



(53) 



+( _ 1)Pn+m [ /n+m ( A )p [g n+m (X)} 



T^n+m — 1 / \ \ 

Li r /\m^2 TT a «+ml A ~ A j 



n 



fe(A _ A (n+m-l) ) > 

while the corresponding Bethe ansatz equations become 

A«)6(A«-A«)-a 2 (Af-A«; 



r/i(\ m )" 


Li 






U(Af'). 




U(A, (1) ). 





^ « 2 (a! 1} 



.^(AW-A^KA^-A^n 6(A f)-A«) 



(54) 
(55) 



15 



\ /*(a? } ) " 


Li 




i 2 






^ + i(Af). 





n 



(fc-l)\ m fc 



(*) 



\ a t (Af» - Af" 1 ') 0t (Af - Af») MA!" - Af) J: 



for 



n + m — 2 



n+m— 1 



(A 



(n+m— 1) \ 



p / X (n+m-l) > 
Jn+mvA y 

( 3'n+m( Aj 



Qn+m—l 



(A; 



(n+m— 1) \ 



/ X (n+m-l)x 



mn+m-2 

n 



6(A 



6(A 



(n+m— 1) 



(fc+1) 



A 



Al fc) ) 



(n+m— 2) \ 



(56) 



n 

3=1 



A 



(n+m— 1) \ 
3 ' 



(n+m— 1) 



3=1 

-A 



G n+m-l(A 



(n+m— 1) 



A 



(n+m-2)> 



(n+m— 1) \ 



/,(n+m-l) 
"n+m-l^j 



A 



(n+m— 1) \ 



6(A 



(n+m— 1) 



A 



(n+m— 1) 



) 



(57) 



Up to this point our analysis has been carried out for arbitrary dressing functions /i(A) and 
<7i(A). In fact, the main results (|53H57|) are expected to be valid for any representations such 
that reference states with the properties (|25H26|) could be exhibited. We note that our results 
for the Bethe ansatz equations (|54H57|) are equivalent to the analyticity of the eigenvalues (J55|) 
as function of the rapidities {A^}, • • • , {A^ ri+m ~ 1 ' ) }. The structure of Eqs. (|54Ho"7|) turns out 
also to be in accordance with that generally postulated in the realm of the analytical Bethe 
ansatz approach [23] for the Sl(n\m) symmetry and its quantum affine extension |2*6l l27j. The 
explicit construction of the transfer matrix eigenvectors as presented here is, however, beyond 
the scope of the analytical Bethe ansatz. 

We now turn our attention to the specific case of the vector and dual representations of 
U q [Sl(n\m)]. In this situation further simplifications can still be implemented in Eqs. (j53l - 
I57j) to bring them to assume a more traditional symmetrical form. This is implemented by 
performing the shift Xj — ► Xj — 5^ where 5® = %\ Ylk=i(~ l) Pfc anc ^ 9 = e ~ K ' ■ Substituting 



the explicit expressions for the weights in Eqs. 



and after some long but straightforward 
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simplifications we find that 
A(A;{A; i) } ) ... ) {Af +m - 1) 



}) = (-ir ma)] Li [b(\+p)ri[ 

i=i 



sinh 


[x? - 


- A - 6<U + i-y( 




sinh 


Aj 1} - A - 





n+m—1 



+ E (-ir[KA)] Li RA+ P )] 



^2 



fc=2 

><n 



sinh 




- 1 ) + i 7 (_i)p fc J 


sinh 


A 


- Af " 1} + tf*" 1 )] 



sinh 


[(Af - 


-A - <SW +i7( 


_1)p* 


sinh 


A« - A - 





^n+m- 1 



+(-l)^[6(A)] L Ma„ +m (A + p)] L2 J] 



n 



sinh [A - +m_1) + !) + 7 (-l)P" 



(58) 



sinh 



and now the non-linear equations for the variables {Aj?} are 



A- A 



(n+m—1) 



_|_ £(n+m-l) 



sinh 



Af } +a(-l) 



pi 



Li 



sinh 



A! 1} -a(-i) 



pi 



m i sinh 
= i sinh 



*f-A« 



A, w - A« - i 7 (-l)P. 



n 



sinh 


[Af 


-Af + 


sinh 


Af 


-Af 



n 



n 



sinh 


Af 






sinh 


Af 


-Af 





sinh 


Af - 


- Af +!7(-l)»+>] 


sinh 


Af 


-Af -i 7 (-l)»] 



k = 2, . . . ,n + m — 2 



sinh 
j=i sinh 



A f+ 1 )_ A f) + ,2(_l)P fc+ i] 



Af +1 >-Af>-ii(-l)^ +1 ]' 



(59) 
(60) 



sinh 


A ( n+m _ 1)+ ^ ( 


_^Pn + m 




[n — m — 2) 


sinh 


A; n+m - i} -^( 


_]_)Pn + m 


+ P-q 


n — m — 2) 



sinh 



* (n+m— 1) \ (n+m— 2) 



X f+m-l) _ x (n+m-2) + ^(_l)P„ +m - 1 



J= i sinh 



•.{n+m— 1) x (n+m— 1) 
Aj — 



+ if(-l)^ 



A i (n+m_1) - / \( n + m - 1 ) _ ^ 7 (_i) P „+ m -i 



L 2 



(61) 
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Note that at the special value p = ry (n^™ _ ij the Bethe ansatz equations (|59H61|) becomes 
a symmetrical system of non-linear equations. We close this section mentioning that the 
eigenvalues £^ of the spin chain H^' discussed in section El are 



"'i 



2sinh[i7(-l) Pl 



i cosh 



2A 



(i) 



where Ajj. satisfy the Bethe ansatz equations 



I- £(0 

cosh [ij(-l) pl ] 

and the additive constants are 



(62) 



C (1) = (L — 1) {coth [ry(— l) Pl ] — coth [i"f]} + {coth [p + ij] — coth [p]} 



(63) 



C (2) = - {coth [i 7 (-l) Pl ] - coth [z 7 ]} + - {coth [p + ij] - coth [p]} . (64) 



5 Conclusions 

In this paper we have studied a trigonometric integrable vertex model that combines the vector 
and dual representations of U q [Sl{n\m)]. The transfer matrix eigenvalue problem was solved 
by means of the algebraic Bethe ansatz method for all possible gradings. From our 

results for the transfer matrix eigenvalues and Bethe ansatz equations one can in principle 
derive the thermodynamic free-energy [2B]- The fact that we know the solution for the many 
possibilities of grading choices could be of great help, since we can choose the particular one 
whose Bethe ansatz roots topology are the less complicated as possible. In particular, it would 
be interesting to determine the classes of universality governing the criticality of the gapless 
regimes. Finding the dependence of the low-lying critical exponents with n, m and 7 may be of 
utility to make connections with supersymmetric Wess-Zumino-Witten conformal field theory 
[T7] . We hope to return to this problem in a future publication. 



Acknowledgements 

The author G.A.P. Ribeiro thanks FAPESP (Fundagao de Amparo a Pesquisa do Estado de 
Sao Paulo) for financial support and L. Amico for the hospitality of the DMFCI, Universita 



18 



di Catania, Italy. The work of M.J. Martins has been support by the Brazilian Research 
Council-CNPq and FAPESP. 



Appendix A: Hamiltonian coefficients 

In this appendix we define the Hamiltonian coefficient in terms of the Boltzmann weights 
defined in section |21 







[Oa(p)ci(p)g- Ma - 1 +Q l {p)ap(j>)( af }] {-l) pp , if a < (3 

[a a (p)t 2 (p)q 2 ^-V + c 2 (p)Mp)l/t Pa \ (-l)W, if a > 

n+m 

Si(p)ik(p)(-lY''**>q- M " 1 Ccn if a <r,l< [3 



7=1 

n+m 



X)SlW^W(- 1 ) P ' ,4 * r 9 2(,A1 " 1) C«r, 



7=1 
n+m 



7=1 
n+m 



V 7=1 



if a < 7; 7 > (3 
if a > 7; 7 < /3 
if a > 7; 7 > (3 



(a.i; 



(i) 



tt {2) 



a,p 



a a (p)c 1 (0)a a (p) if a > f3 
a a (p)c 2 (0)d a (p) ifa<(3 

« a (p)«a(0)ci(p)(-l) Pa+p ' 3 g- 25Q ' 1 ifa<p 
a (p)a a {0)c 2 (p){-iy«+^q 2 ^-V if a >(3 

^(j3)ap(0)ap(jj)Cap if a < (3 
c 2 (p)dp(0)ap(p)l/Cp a ifa>(3 



(A.2) 



(A.3) 



(A.4) 
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q a (p)c 1 (0)c 1 (p)(-l 


|P/3g-2<5 a ,l 




o a (p)c 1 (0)c 2 (p)(-l 


1^/3^2(5,3,1-1) 




a a (p)c 2 (0)c 1 (p)(-l 


|P/3g-2<5 a ,l 




o a (p)c 2 (0)c 2 (p)(-l 


|P/3g 2 (5^,i-l) 




c 1 (p)a 7 (0)ci(p)(-i; 






c 2 (p)a 7 (0)ci(p)(-i; 






c 1 (p)a 7 (0)c 2 (p)(-i; 






c 2 (p)d 7 (0)c 2 (p)(-i; 




ya 


c 1 (p)c 2 (0)a j9 (p)(-l 






c 1 (p)ci(0)a /3 (p)(-l 






c 2 (p)c 2 (0)a /3 (p)(-l 






< c 2 (p)ci(0)a«(p)(-l 

— z Vr / J- \ / p / \ 

n-\-m 




Sew 


v=\ 

n+m 

^c 1 (p)c 2 (0)c 2 (p) 


r _iyPf3+Pvq2(Sf3, 


L-1)A 

Sai/ 


u=l 
n+m 

c 1 (p)c l (0)c 1 (p) 


' —lyp+P" q~ 2S " 




u=l 

n+m 

^c 1 (p)c 1 (0)c 2 (p) 




L-1)A 


u=l 

n+m 

c 2 (p)c 2 (0)ci(p) 


' _\yi3+Pvq~ 2 5 v 




n+m 

c 2 (p)c 2 (0)c 2 (p) 


'—iypp+Pvq 2 ( 5 0, 




I/=l 

n+m 

Yl c 2 (p)ci(0)ci(p) 


' —iy0+P"q- 2S « 




n+m 

5^ c 2 (p)ci(0)c 2 (p) 


'_X)P/3+P"g2(<5 (3i 


l - 1} vc 



if a < (3 

if a > (3 

if a < (3 

if a > (3 

if o; < 7 

if a > 7 

if a < 7 

if a > 7 

if a < /3 

if a < (3 

if a > (3 

if a > f3 

if a < v 

if a < v 

if a < v 

if a < v 

if a > v 

if a > v 

if a > v 

if a > v 



a < 7 
a < 7 
a > 7 
a > 7 

0> 7 
/3> 7 
(3< 1 
(3< 1 
(3> 1 
(3< 1 
(3> 1 
(3< 1 

(3 > v\ 7 < v 
f3 < h>; 7 < v 
f3 > h>] 7 > v 
(3 < v; 7 > v 
f3 > h>; 7 < v 
(3 < v; 7 < v 
f3 > h>; 7 > v 
f3 < h>] 7 > v 



(A.5) 
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^(1) 



,(2) 



,(3) 



,(4) 



6(p)ci(0)6(p) if a < (3; a ^ 7 
b(p)c 2 (0)b(p) if a > f3; a ^ 7 

a a (p)6(0)6(p)(-l)^5 a)7 ifa^/J 
fe(p)&(0)^(p)(-l)^^, 7 

6(p)6(0)6(p)(-l) PaP " if a ^ (3; a ^ 7; P + 7 

6(p)6(0)c 1 (p)(-l)P^+P' 3 g- 2 ^. 1 if a ^7; /3 > 7 
&(p)&(0)c 2 (p)(-l) p ^ +p V (5/M_1) if a ^ 7; /3 < 7 

c 1 (p)6(0)6(p)(-l)^w+^Ca 7 if a < 7; /3 7^ 7 
c 2 (p)6(0)6(p)(-l)^+^l/C 7Q if a > 7; /? ^ 7 



(A.6) 



(A.7) 



(A.8) 



(A.9) 
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